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Exact diagonalizations with a realistic interaction show that configurations with four neutrons 
in a major shell and four protons in another — or the same — major shell, behave systematically as 
backbending rotors. The dominance of the q ■ q component of the interaction is explained by an 
approximate form of SU3 symmetry. It is suggested that these configurations are associated with 
the onset of rotational motion in medium and heavy nuclei. 
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, The SU3 model of Elliott |1J] provides a microscopic description of rotors that exhibit spectra in J(J +1). For 
Q\ • sufficiently low J, or sufficiently large representations they became perfect in the sense of having a constant intrinsic 
^ — 1 , quadrupole moment Qq = Qq(J), where 

qq ' as postulated in the strong coupling limit of the unified model of Bohr and Mottelson . 

Since the quadrupole force that appears in the SU3 Casimir operator is also an important part of the nuclear 
interaction ||Q , we expect it to play a determinant role in the onset of rotational motion in real nuclei - the problem 
we want to address. A direct approach would demand, in general, diagonalizations in spaces of two major shells in 
neutrons and protons as first proposed by Kumar and Baranger H. Dimensionalities are then of order 10 40 , exceeding 
. by far what is possible at present (10 7 ) ||. 
I/"") ' Therefore, it is necesssary to develop a computational strategy, and our starting point will consist in learning 
as much as we can from situations in which neutrons and protons are independently restricted to a single major 
shell, that can be the same close to N=Z. The exact results in I will show that rotational features -including the 
systematic appearance of backbending- are determined by the interplay of the quadrupole force with the central field, 
in the subspace of a major shell spanned by the sequence of Aj = 2 orbits that comes lowest under the spin-orbit 
splitting. This state of affairs is related to the existence of an approximate symmetry (quasi-SU3) , introduced in II. 
' The relevance of these results to rotors in medium and heavy nuclei will be explained in III. 
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NOTATIONS . v = neutrons, n = protons, C lm = ^/4tt/(2Z + l)Y lm , q = q 2m = r 2 C 2m . 

. p is the principal quantum number, r p is the generic label for all orbits in the p-th oscillator shell except 

the largest (i.e. j = j nmx =p + 1/2). 

1_h ' We use I for j = I + 1/2 orbits in the sense h = /in/2, g = <?g/2j V — P3/2 etc., except in the following 

CONVENTION: pfh means the full p = 5 shell, i.e. £>i/ 2 ,f>3/2> fs/2,, St/2, ^9/2, h n / 2 , while hfp = 
^■11/2/7/2^3/2! an d similarly for other shells. 

I. Exact Results. Although a space of a full major shell, with very specific single-particle spacings, is necessary 
to ensure strict SU3 symmetry, we know of several examples where the ds or fp subspaces produce rotor-like spectra 
in the presence of spin-orbit splittings: (ds) 4 describes 20 Ne quite well M and (ds)%' A (f p) 2 , configurations explain the 
onset of deformation in 31 Na and 32 Mg S. Furthermore 48 Cr provides the first example of a backbending band in 
N=Z nuclei. The experimental spectrum M is almost perfectly reproduced by a full (pf) s shell model calculation, 
with strong indications that the (fp) s space is sufficient to explain the quadrupole coherence The situation has 
a double interest. As we shall explain in III, configurations that consist of 4 protons in a major shell and 4 neutrons 
in another (the same in NwZ) play a key role in the onset of rotational motion in heavier nuclei, and the restriction 
to the Aj = 2 spaces makes the diagonalizations possible, as illustrated by the four cases we are going to treat (in 
parenthesis the corresponding m-scheme dimensionalities): 
(f P ) s T = ; (2 x 10 4 ), (fp)U9ds)i ; (1.1 x 10 5 ), 

(gds) 8 T = ; (6 x 10 5 ), (gds)i(hfp)i ; (1.9 x 10 6 ), 
against 

(pffT = 0; (2xl0 6 ), (Pf)i(sdg)t ; (10 7 ), 

(sdgfT = ; (5 x 10 7 ), (adg)*(pfh)i ; (1.9 x 10 8 ). 
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We shall compare the results obtained with the KLS interaction jl2[ and with a pure quadrupole forces using 
Uu) = 9 MeV with a uniform single particle spacing e — 1 MeV, corresponding to the standard — (3d ■ s splitting 
(13 w 20A- 2 / 3 MeV and Tilu w 40 A' 1 / 3 ) jb|. 

It has been shown in [[| that for one shell the quadrupole component of a general realistic interaction has the form 
— e2<7 P ■ <7 P , where e-i goes as A -1 / 3 and, q p — q p /J\f p 2 ' 1 is the quadrupole operator q p in shell p, normalized by Mp - 
the square root of the sum of the squares of the matrix elements of q p - which goes as (p + 3/2) 2 . For two contiguous 
shells the force is —e2(q p + Qp+i) ' (%> + Qp+i)i with the same e2 coupling, and it differs markedly from the traditional 
x(q P + q P +i) ■ {q P + q P +i), with x = 0(A- 5 / 3 ) f§. 

Fig.l shows the yrast bands in the four spaces. Rotational behaviour is fair to excellent at low J. As expected 
from the normalization property of the realistic quadrupole force the moments of inertia in the rotational region go 
as (p + 3/2) 2 (p' + 3/2) 2 , i.e. if we multiply all the E 1 values by this factor the lines become parallel. The Qo values 
are constant to within 5% up to a critical J value at which the bands backbend. 

Since all the spaces behave in the same way we specialize to (gds) 8 in what follows. Fig. 2 shows the results of 
diagonalizing e<i q p ■ q p (p = 4). At e-i = 9.6 the s splittings are overwhelmed and we have a nearly perfect rotor. The 
value of Qo stays practically constant up to J = 16 — 18 and then decreases slowly. At = 4.8, 3.2 and 2.4 the 
rotational behaviour remains very good below J = 14. Then there is a break and the upper values are again aligned. 
At e2 = 3.2 the overlap of each state with the one obtained with the full KLS interaction is always better than (0.95) 2 , 
which suggests that 

< h\H\h >,/«< q\U\q >,/, 

where \h > and \q > are the eigenstates of the full Hamiltonian H and the quadrupole force [e^ — 3.2) respectively. 
Fig. 3 shows that this is the case indeed. It means that the observed backbending pattern is obtained by doing first 
order perturbation theory on \q >: the spectrum changes but not the structure (i.e. the wavefunctions). A similar 
situation is found when comparing the full (pf) s calculation with a renormalized interaction and e = 2 (fig. 10 of ref. 
JTT||) and the (fp) s result in fig. 1: the backbend occurs at the same J and the Qo values are very close in spite of a 
much larger slope (i.e., smaller moment of inertia) in the bigger calculation. Here again perturbation theory should 
operate well. 

To gain some insight into the backbending phenomenon we examine the evolution of the wavefunctions and 
quadrupole moments. In fig. 4 we find that for e2=9.6 the percentage of the g s configuration in the full eigen- 
state is very small and nicely correlated with the Qo values. This is what we expect from a good rotor, for which the 
amplitudes of any configuration (not only g s ) must be J-independent (since all states must be projections of the same 
intrinsic state). For the KLS results and their e2=3.2 counterparts Qo(J) decreases abruptly above J = 14, while the 
g & configuration increases its amplitude and becomes dominant in the region where Qo(^) reaches a plateau. It is 
clear that at the backbend the notion of intrinsic state loses - or changes - its meaning, and the idea of a band crossing 
suggested by fig. 2 becomes questionable. Work remains to be done to understand the connection of our results with 
the mean field ones fl4|| . 

II. Quasi- SU3. That the build up of quadrupole coherence needs only the lower Aj = 2 sequence of the full 
shell can be understood by examining table 1 where we list the matrix elements of q 20 — r 2 C 20 — \{Sz 2 — r 2 ), in 
jj and LS coupling. It is seen that the Aj = 1 matrix elements are small, both for to small (prolate shapes) and m 
large (oblate). If we simply neglect them, diagonalizing the Aj = 2 matrix in jj scheme is very much equivalent to 
diagonalizing the exact q 20 operator in LS scheme. This amounts to saying that the sequence j = 1/2, 5/2, 9/2 • • • 
(or j = 3/2, 7/2, 11/2 • • •) must behave very much as an I — 0, 2, 4 • • • (or I = 1, 3, 5 • • •) one. Therefore we introduce 
a new operator (the 'quasi' q 2a ), defined in the Aj=2 space via the following replacements in the LS matrix elements 
of q 20 : I — > j , p — > p + 1/2, to — > rn + 1/2 and — m — > — m — 1/2: (to > 0). 

In fig. 5 we draw to the left the spectrum of the full q 20 operator (in fact 2q 20 ), i.e., the SU3 Nilsson orbits. The 
band-heads come at 2(p + 3/4 — 3/2|to|). To the right we have plotted the spectrum of the 'quasi' 2q 20 operator. Now 
the band-heads are at 2(p + 1/2 — 3/2|m|), that is, the exact LS values, except for m = ±1/2, where the one to one 
correspondence between the "quasi" g 20 and the exact q 20 in LS scheme breaks down. The corresponding "quasi-SU3" 
symmetry cannot be exact because of this mismatch. Notice however that the mismatch is very small (< 1%). The 
spectrum of the true q 20 operator in the Aj = 2 space is extremely close to the one in fig. 5, and it is clear that the 
amount of quadrupole coherence obtained by filling the m (or K) — 1/2 and 3/2 orbits is almost as large as for the 
SU3 orbitals. For the eight particle blocks we are interested in, the intrinsic Qo would be 

Qo = 8M^ - 1) + e v (jp v - 1)] (2) 

For p n = 5, p u = 6 and effective charges e n = 1.5, e„ = 0.5, we have Qo = 68 (in dimensionless oscillator coordinates, 
i.e., r —* r/b with b w A~ x / 6 /to), which is related to the E2 transition probability from the ground state by 
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B(E2) |= 10- 5 A 2 / 3 Qg = i.4 e 2 fe 2 for A = 166 We shall 

see soon the relevance of this number. The Qo values 
obtained in I with the q ■ q interaction at e-i = 9.6 saturate the value predicted by eq.2 within 2% while at e2 = 2.4 
we still have 80% of this limit. 
III. The onset of rotation. 

Fig. 6 proposes a schematic view of the single particle order expected around Z = 50, N — 82 closure. Mutatis- 
mutandis, the same scheme applies to the Z = 28, N — 50 and Z = 82, N = 126 closures. For the lower shells we 
have a conventional sequence, and for the upper (empty) ones we have assumed a spin orbit splitting, which may be 
naive, but it is correct in the light nuclei and consistent with the (scarce) available data in the heavier ones . 

Nilsson diagrams |l5| , |l6| ] predict that when nuclei acquire a stable deformation, two orbits K = 1/2 and 3/2 
- originating in the upper shells of fig. 5 - become occupied. This result is common to different regions and to 
different calculations and provides a strong clue that we turn upside down when translating it into spherical language 
as: Rotational motion sets in when 8 particles are promoted to the (pfh)^(sdgi)f, configuration, suggesting that the 
calculations we have presented may be the first step in implementing a shell model view of real rotors. The BE2 
estimate at the end of the previous section is a factor four smaller than the largest observed values in the region 
fl7|| . It means that Qo is only a factor two too small. Therefore the lower orbits should supply the missing half, 
which is very plausible because their potential quadrupole coherence is high. The reason is that the r p groups in 
fig. 6 are pseudo-oscillator shells with p' = p — 1 and the pseudo-SU3 symmetry of Arima, Draayer, Harvey and 
Hccht (and hence the left part of fig. 5) would provide the relevant coupling scheme (ref. |l^] contains a recent survey 
of pseudo-SU3). It should be noted that the quadrupole coherence in the upper(u) and lower(l) spaces is mutually 
reinforced through coupling terms of the form —eiq u ■ <7z- 

Conclusions. The zeroth order approximation we are proposing for rotational nuclei is very similar to the weak 
coupling explanation of the famous 4p-4h low lying states in 16 O and 40 Ca. The differences are that in the heavier nuclei 
they become ground states, as in 80 Zr, and most probably 8p-8h excitations are necessary to ensure the observed 
quadrupole coherence. The fact that the calculations naturally produce backbending in the upper configurations 
(which are the driving ones, as in light nuclei) indicates that we are closer to real rotors than to idealized constructions. 
The tendency of quadrupole forces of Elliott type to produce clustering in the excited states [ 
lead to significant differences of interpretation between the spherical and deformed formulations for large quadrupole 
moments. 

This work is supported in part by the IN2P3 (France)-CICYT (Spain) agreements an by grant DGICYT PB93-263 
(Spain). 
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TABLE I. The matrix elements of r 2 and C20 in jj and LS coupling 



<pl\r 2 \pl >=p + 3/2 



< jm\Ci\jm >- 



j(j + 1) - 3m 2 



< jm\C2\j + lm >= 



2j(2j + 2) 

3m[(j + l) 2 -m 2 } 1 ' 2 



(2j + 4)(2j + 2)(2j) 



< jm\C 2 \j + 2m> 



_ 3 ( [(j + 2) 2 - 
2 1 (2j + 



m 2 ][(j + l) 2 -mV ' " 



2) 2 (2j + 4) 2 



<pl\r 2 \pl + 2 >= -[(p-/)(p + « + 3)] 1/2 
J(Z + 1) - 3m 2 



< Zm|C 2 |Zm >= 



(2Z + 3)(2i — 1) 



< lm\C 2 \l + 2m >= 



3 j [( ? + 2) 2 -m 2 ][(l + l) 2 -m 
(2Z + 5)(2Z + 3) 2 (2£ + l) 



21 1 V2 



4 




FIG. 1. Yrast transition energies E n = E(J + 2) — E(J) for different configurations, KLS interaction. 




FIG. 2. Yrast transition = E(J + 2) — E(J) for the (gds) 8 configuration with an ~e-zq ■ q force. 
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FIG. 3. < h\H\h >= (gds) 8 in fig. 3; < q\qq\q >= 3.2 in fig. 4 compared with < q\H\q >. See text. 




FIG. 4. Qo{J)/Qo{2) (full lines) and g$ =< g 8 \(gds) 8 > 2 (dashed lines). Wavefunctions calculated with — e^q ■ q (crosses 
= e2 — 9.6, squares = e2 = 3.2) and KLS (diamonds). 
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FIG. 5. Nilsson orbits for SU3 (k = 2p) and quasi-SU3 (k = 2p - 1/2). 
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FIG. 6. Schematic single particle spectrum above 132 Sn. 
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